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ABSTRACT 



Context. We present 12.8 jim images of the core of NGC 1068, the archetype Seyfert type II galaxy, obtained during first operations of the 
BURST mode of the VLT/VISIR (Imager and Spectrometer in the InfraRed at the Very Large Telescope). 

Aims. We trace structures under the diffraction limit of one UT (Unit Telescope at the VLT) and we investigate the link between dust in the 
vicinity of the central engine of NGC 1068, recently resolved by interferometry with MIDI (Mid-InfrareD Interferometer), and more extended 
structures. This step is mandatory for a multi-scale understanding of the sources of mid-infrared emission in active galactic nuclei (AGN). 
Methods. A speckle processing of VISIR BURST mode images was performed to extract very low spatial-frequency visibilities, first considering 
the full field of VISIR BURST mode images and then limiting it to the mask used for the acquisition of MIDI data. 

Results. Extracted visibilities are reproduced with a multi-component model. We identify two major sources of emission at 12.8 yum: a compact 
one < 85 mas, directly associated with the dusty torus, and an elliptical one of size (< 140) mas x 1 187 mas and PA. ~ -4° (from north to east), 
which gives a new description of the NS elongation of the nucleus. This is consistent with previous deconvolution processes. 
The combination with MIDI data reveals the close environment of the dusty torus to contribute to ~ 83% of the MIR flux seen by MIDI. 
Conclusions. This strong contribution has to be considered in modeling long baseline interferometric data. It must be related to the NS 
elongated component which is thought to originate from individually unresolved dusty clouds and is located inside the ionization cone where it 
is photoevaporating and radiatively accelerated. Low temperatures of the dusty torus are not challenged, emphasizing the scenarios of clumpy 
torus. 
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1. Introduction 

NGC 1068, a bright and near archetype Seyfert type II galaxy 
is a key target in testing the unified scheme of active galac- 
tic nuclei (AGN). This stipulates that Seyfert II and Seyfert 
I are the same objects, depending upon the inclination of th e 
nucleus along the line of sight (lAntonucci & Millen Il985l) . 
The dusty torus is therefore an important element of this 
paradigm. The most r ecent theoretical models a re given by 



upper size limits of 50 mas in the L' and M' bands were set 



(iGratadour et al.ll2006h . These data also exhibit an ~ 800 mas 



Nenkova et al] (|2002|) 



(120051) . iHonig et al 



.Schartmann et al . 

(l2006h . lElitzur & Shlosmanl (l2006b . among others 

Multi-wavelength observations of t he nuc leus of NGC 1068 
have been reviewed bv lGalliano et al.l (l2003b . Recently, decon- 



volution of near-IR adaptive optics (AO) images of NGC 1068 
shows an unresolved core of some 30 mas along the NS axis 
and < 15 mas along the EW axis (in the Kj band). Additionally, 



* Based on commissioning time observations collected at the 
ESO/Paranal MELIPAL telescope. 



extended, peculiar S-shaped distribution aligned with the ra- 
dio jet. Near-IR speckle interferometry has revealed a com- 
pon ent of FWHM ~ 1 8 mas x 39 mas oriented at PA ~ - 
16° (feigelt et al."2004'). Combined with the VLT/VINCI mea- 
surement of (iWittkowski et al..,2004) . data were fitted with two 
components of F WHM < 5 mas and 40 mas. The first MIR im- 
ages reported by iBock et al. I (l2000l) traced the same NS elon- 
gation of the nucleus, and the 12.8 um de convolved images 
from the VLT/VISIR jGalliano et al.ll2005b highlights several 
knots distributed along the S-shape (twis ting from PA ^ -4° 
to PA = 31°). Using mid-IR spectroscopy, iMasoneLaL I (l2006l) 
distinguished two sources of mid-IR emission, a compact core 
within the central 4"x4" and a more extended component 
linked to heated dust inside the ionization cone. Mid-IR in- 
terferometric observations with the VLT/MIDI were necessary 
to resolve the dusty torus of NGC 1068. The first data were 
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obtained in 2003 and presented by iJaffe et alj (120041) . Two in- 
dependent analyses have revealed a compact structure of dust 
composed of amorphous silicates, with size < 82 m as and rel 



atively low temperatures (on the order of 350 K; iJaffe et al 
2004l:|Poncelet et al.ll2006l) . 



Configurations of Unit Telescopes (UTs) at the VLTI do 
not provide baselines shorter than 30 m, although low spatial- 
frequency visibility points would yield tight constraints for 
modeling MIDI data. One way to at least partially solve this 
problem is observe with a single UT to access baselines be- 
tween and 8.2 m (the pupil of one UT) using an interferomet- 
ric processing of images. 

NGC 1068 was observed at 12.8 fim in January 2005 using 
the BURST mode (or quasi-Speckle mode) of VISIR. The big 
advantage of this is the shorter acquisition time of 16 millisec- 
onds, far shorter than the typical atmospheric turbulence varia- 
tions of ~ 300 ms at 10 jim. The aim of the study is therefore 
to present the BURST mode of VISIR to show its advantages 
and its use for the speckle analysis. This allows us to derive the 
mid-IR low spatial-frequency visibilities between and 8 m of 
baseline for the first time. Next we compare these to MIDI data 
points in order to establish the link between the source of mid- 
IR emission on different scales. Section |2] gives details on the 
BURST mode of VISIR and presents the principle of how visi- 
bilities are extracted from VISIR images. Results of the model- 
ing applied to the VISIR visibilities are shown and discussed in 
Sect. E] Section Hlpresents the Hnk between VISIR and MIDI 
observations. Section|5]gives the general conclusions drawn by 
this study. 

2. Observations and data processing 

NGC 1068 was observed on the night of January 25, 2005 us- 
ing the BURST mode of VISIR, which is installed on the VLT 
(P aranal, Chili). The actual observation mode used is described 
in 



Doucet et al 



2006). Total exposure time on the object was 
208 seconds using the 12.8 yum Ne II filter (A^ = 0.21 /im). The 
exposure time was 16 ms for individual images. The airmass 
was ~ 1.20 at the time of observations. HD 11383, a K5III spec- 
tral type star, was observed as a calibrator. To correct for the 
turbulence with offline processing, data are stored every 1000 
elementary images for one nodding position, using a chopping 
frequency of 0.25 Hz in the NS direction. After classical data 
reduction in the mid-IR range, a cube of 500 images, both 
chopped and nodded (4 beams/image), is obtained. Turbulence 
has the effect of moving each of the 4 sources independently so 
they have to be individually extracted from the cube. Finally 
images of similar quarter are shifted and added, producing 
some 13000 and 2200 elementary images of NGC 1068 and 
HD 11383, respectively. To obtain a better signal-to-noise ra- 
tio - and according to the good correlation between elemen- 
tary images - they are added in blocks of 5 (the sums may 
be considered as instantaneous maps of the source intensity). 
Individual and resulting 12.8 ;um VISIR BURST mode images 
ofNGC 1068 andofHD 1 1383 may be found in Fig. [U 

The BURST mode is important because it uses such a short 
exposure time that it effectively freezes atmospheric turbu- 
lence. Therefore it enabled us to apply the speckle technique 



first introduced by iLabeyrie (Il970l) as an analytical tool. This 
aims to reach the diffraction limit of a large telescope (i.e. 
~ 320 mas at 12.8 mas), after incoherent integration of expo- 
sures. 

The seeing was excellent and very stable during observa- 
tions of NGC 1068 and images are diffraction-limited. Thus, 
the monochromatic transfer function f ,i(M, v) is estimated well 
by observations of the calibrator star Unbiased visibilities for 
NGC 1068 and for the calibrator star are computed from the 
power spectra of the images: 



|Vcahb(", V)P = <|/cahb {U,V)\\ ' (l^cahb (u,v)\\ (1) 

|VnGC1068 (", v)|^ = <|/nGCI068 (m, v)|^), - (lA^NGCioes (", v)|^), 

where /ngcio68, /caiib, A^ngcio68, and ./Vcaiib are the Fourier 
transforms of the intensity distribution on individual VISIR 
BURST mode images of NGC 1068, of the calibrator star, and 
of the background photon noise on the related images, respec- 
tively. The final calibrated visibility estimate is 



w 



NGC 1068 



(",v)p 



IVcalib (", V)P 



(2) 



In our study, the most important interest of BURST mode 
observations compared to long-time exposures is the ability 
to calibrate the source spatial spectrum. We considered error 
propagation to determine error bars for the uncalibrated vis- 
ibilities. A Monte-Carlo computation is used to estimate er- 
ror bars for the calibrated visibility. It consists in generating 
random numbers with a Gaussian distribution for IVNGCioesP 
and |y caiibP (as the histograms of the individual exposure esti- 
mates are Gaussians) taking their standard deviations into ac- 
count. The standard deviation of the simulated visibility distri- 
bution is the 1 cr error on \V{u, v)p. Depending on spatial fre- 
quencies, error bars on visibilities range between 1 and 10 %. 
Two-dimension (2D) maps of the visibilities and associated er- 
ror bars are obtained. Because the intensity distribution of the 
source is real, this 2D visibility map is centro-symmetric with 
respect to the zero spatial frequency (see left of Fig. |2]i. Cuts 
of the map along several orientations are presented in Fig. [3] 
Orientations of the spatial frequency plane correspond to the 
orientation on the sky plane (conventions used being 0° is re- 
lated to the NS axis and 90° to the EW axis in images). 

In the following, the visibilities are compared to simple 
geometrical models as is often the case in optical-infrared in- 
terferometry. A model of the spatial brightness distribution of 
the source 0{a, 6) is used to extrapolate the high-spatial fre- 
quency information from the measured low spatial frequencies. 
Thus, thanks to the a priori information injected through this 
model, structures smaller than the diffraction limit of the instru- 
ment can be estimated. This super-resolution technique works 
as long as the signal-to-noise ratio on visibilities is high enough 
to be able to detect a visibility drop with respect to 1 . In clas- 
sical imaging, an equivalent technique consists in detecting the 
width difference between the point spread function (psf) and 
the source intensity peak. This requires a high signal-to-noise 
ratio, however, and an excellent stability of the psf. 



Poncelet et al.: An original interferometric study of NGC 1068 with VISIR BURST mode images 



3 




Fig. 1. 12.8 fjm images obtained with the BURST mode of VISIR. For comparison, normalized images are displayed with the 
same color table. From left to right: a) Sum of the 2200 individual exposures on HD 1 1383 (reference star). The good signal-to- 
noise ratio has to be compared to individual exposures. Here we clearly see the Airy disk (~ 10.5 pixels large centered on the 
brightest pixel), b) Sum of the 13000 exposures obtained on NGC 1068. The full structure appears clearly: the NS elongation 
from the unresolved compact core, the asymmetric emission and a knot farther away toward the NE (north is up and east is 
left). The appearant EW elongation is due to the diffraction ring at ~ 393 mas from the core, c) Example of an individual 16 ms 
exposure of HD 1 1383. The Airy disk faintly appears, attesting to the good seing conditions during the night of observation, d) 
An individual exposure of NGC 1068, already showing the NS elongation of the nucleus. 



3. First study of VISIR visibilities 

In this section, visibilities are compared to uniform disk mod- 
els that we used as rulers for measuring the typical sizes of 
structures in the object as a function of azimuth. The complex- 
ity of models increases to account for the whole evolution of 
the visibility and especially for the departure from a circular 
symmetry. 

3.1. Modeling 

3.1 .1 . Azimuth-dependent modeling 

The 2D visibility map derived from VISIR images (see left of 
Fig|2]i shows an elongated shape in a direction close to 90°. 
Similarly, cuts of this map along different orientations (see 
Fig. [3]) highlight the steep fall in visibility along the NS axis 
(see plots corresponding to 0° in Fig.O in comparison with the 
perpendicular axis (see plots corresponding to 90°). The higher 
the visibility, the smaller the observed source. The source is 
then well-resolved along the NS axis and smaller as we look 
towards the perpendicular axis. A minimum of two uniform 
disks in the model is necessary to fit the different cuts of vis- 
ibilities. For each orientation (marked by 6), the free parame- 
ters associated with the model are the sizes of the two disks 
(01 (0) and 02(0)) and the flux ratio between them {7]{9)). Fits 
are presented in Fig. [3] Along each direction, optimum param- 
eters were evaluated with a x~ minimization and error bars on 
parameters were obtained by varying the ■ by 1 . The global 
value of the reduced (i.e. corresponding to all fits) is given 
by 

2 _ 1 x^\Vf-Mi{0i{9), 02(0), 7]{9) )f 



where (A^ - Pi) is the number of degrees of freedom, stands 
for the number of visibility points (A^ = 91), and Pi is the 



number of free parameters. There are 3 parameters per orienta- 
tion and a total of 12 orientations (Pi = 36), V, and cr, are re- 
spectively visibilities and error bars measured with the speckle 
analysis, M, are the modeled visibility values. The value of the 
global reduced is 24.9. The size of the inner component is 
< 300 mas, below the diffraction limit of the UT. The second 
component is well-resolved (~ 1600 mas) and more or less cir- 
cularly symmetric. For each component, sizes and error bars 
related to each orientation are reported in the image and linked 
for the extended one (left panel in Fig.|4]i. Nevertheless, it is dif- 
ficult to figure out the real shape of the components since part 
of the asymmetry is due to the azimuth dependence of the flux 
ratio between the two components. It ranges between ~ 0.29 
along the EW axis and ~ 4.76 along RA. = 14° from N to E. 
A more realistic description of MIR sources in the nucleus of 
NGC 1068 is given by the global modeling explained below. 

3.1.2. Global modeling 

Using our previous results, we attempted to produce a global 
fit of all orientations at the same time using a more complex 
model. The global fit of the 2D map of VISIR was based on the 
most reliable visibility points (see left of Fig|2]i. We here con- 
sider the flux ratio between the components to be a constant. 
To account for the asymmetry of the visibility distribution, we 
first made use of a simple model of one circular, plus one ellip- 
tical, uniform disk. This model converges to two different kinds 
of solutions with an equivalent value of the reduced ;\^-=27. 
This value is equivalent to the one obtained with the previous 
model applied to each direction independently. It turns out that 
an additive component was mandatory to have a better fit of 
the visibility map and to better constrain the upper size limit 
of the smallest component. To limit the number of parameters 
and account for the asymmetry of the visibility map, this third 
model consists in two circular uniform disks plus an elliptical 
one. There are seven free parameters: the diameter of the two 
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Fig. 2. Left: visibility map as extracted from BURST mode VISIR images. Only data points used for fits have value, others are 
set to 0. Middle: visibility map from the model of two circular disks and one elliptical uniform one (reduced ;\f^=8. 7). These two 
maps are scaled between and 1 (white points corresponding to a visibility value of 1 and black ones to value of 0). Comparison 
between left and right illustrates the quality of the fit. Right: map of absolute terms of residuals, i.e. substraction of the modeled 
visibiUty map to the measurements one. The color scaling of the residuals map is presented on the right side. 



circular disks (0i and 02); the semi-major and the semi-minor 
axes a and b, and the orientation of the elliptical disk in the 
sky plane (angles are positive from north to east) ; the flux ra- 
tios between the small and intermediate components (771) and 
between the intermediate and the largest one (772). 



2 



1 



N rv2 



N-P2 



Vf - Mi ( 01, 02, a, b, 6, 771, 772 ) 



(4) 



where = 91 and P2 = 7. 

The minimization of the previous expression of the reduced 
yields the following constraints on the parameters; 

- 01 < 85 mas 

- a < 140 mas 

1187 + 22 mas 
-0 = -3.7 + O.8° 

- 02 = 2367 + 107 mas 
-771 =3.0 + 0.1 

-772 = 1.4 ±0.1 

Values for 61 and a are close to zero. We chose to give an 
upper limit for these parameters as is within the 1 cr range of 
uncertainty. The corresponding modeled visibility and residu- 
als maps are presented in the middle and on the right of Fig|2] 
respectively. Fits of visibilities by this model for different ori- 
entations in the (u, v) plane are compared to the ones obtained 
with the first model of two circularly uniform disks in Fig. [3] 
The reduced is equal to 8.7. It is smaller than with the first 
model for which the number of degrees of freedom was much 
larger 

For comparison with previous observations, geometrical 
parameters const rained by this model a re sh own on the decon- 
volved image of lOaUiano et al.l (l2005b and iBock et"ai] (I2OOOI) 
(see the middle and right panels of Fig. HJ. The compact core 
and the ellipse account for the core and the NS elongation (i.e. 
knots NEl and SWl of lGaUiano et al.l2005b . while the 2'.'4 uni- 
form disk accounts (both in size and total brightness) for the 
background of fainter structures at larger scales. The compari- 
son shows the consistency of the different processings applied 
to MIR images obtained at Keck 2 and with VISIR at the VLT. 



3.2. Discussion of tine results 
3.2.1. The MIR core 

The speckle analysis provided a constraint in size, i.e. an up- 
per limit of ~ 85 mas for the inner source of emission in the 
nucleus. This is well below the diffraction limit of the UT at 
12.8 //m (i.e. 320 mas) or the FWHM upper Umits of 290 
and 190 mas (along the NS and EW directions respectively), 
which are measured from deconvolution of the 12.8 fim im- 
age obtained with t he standard acquisition mode of VISIR 
(iGalhano et al.ll2005h . 



This structure is directly associated with the dusty torus 
first resolved by interferometry using MIDI. It is in g ood agree 



ment with sizes produced by studies from both [jaffe et al . 
( |2004 FWHM < 49 mas) and lPoncelet et al.1 ( l2006l diameter 



< 82 mas). The present study does not allow any questions to be 
addressed about the orientation or inclination of this dust struc- 
ture. However, evidence fro m H2O maser disks inclined at -45° 
( Greenhill & Gwimil 1997 ) is though t to indicate the lo cation 
of the inner edge of the dusty torus jKartie et al.lll999b . This 
orientation has subsequentl y been confirmed by the analysis of 
recent MIDI observations (IJaffe et al.ll2006l) . so we adopted it 
in Fig.|5] 

In the NIR, the K-band core has be en found to measur e 
18 mas X 39 mas oriented at PA ~ 16° dWeigek etal.ll2004 . 
This orientation differs from both the MIR and the distribution 
of water masers. Like wise, deconvolved NACO images from 
iGratadour et al. (2006) trace structures ~ 30 mas and < 15 mas 
in the K-band along the NS and EW axes, respectively. This 
K-band structure is thought to be associated to the inner edge 
of the dusty torus, which is heated by radiation from the central 
engine. 

3.2.2. The NS elongated structure 

The (< 140) mas x 1 185 mas elliptical structure is aligned with 
the base of the radio jet in the AGN of NGC 1068 (P.A. ~ -4°). 
The ellipse matches both the NEl and S Wl knots of the de- 
convolved image of iGalliano et al.l (|2005[ see middle panel in 
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Fig. 3. VISIR visibilities derived from BURST mode images at 12.8 fim and associated error bars (symbols) along several direc- 
tions (angles are positive from north to east). One can see the softer fall of the visibility toward the EW axis (90°) than toward the 
NS axis (0°), highlighting the NS elongation of the source. Black solid line corresponds to fits with the two circular uniform disk 
model independently applied for each orientation (global ~ 24.9). Dashed grey line corresponds to the last model considered 
consisting in a two circularly and one elliptical uniform disk. This model is used in order to fit all baselines simultaneously. The 
value of the reduced is here 8.7. It is smaller than the value related to the previous model due to the higher number of freedom 
degrees of the global model (see Sect. 13. 11 1. 
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Fig. 4. Left: comparison of different contributions to MIR emission in the core of NGC 1068. Background is the resulting 12.8 
yum VISIR BURST mode image. Sizes derived from the two uniform disks model applied to each orientation and as described in 
Sect. l3.1.Tl are superposed {red and green components). The black square displays the MIDI mask considered in the second part 
of the study of VISIR visibilities (see Sect.|4|i. Middle: comparison between the structures resulting fro m the global modelin g 
presented in Sect. 13. 1.21 and the deconvolved image obtained from VISIR standard-mode observations ("Galliano et aL 2005 ). 
Right: same compa rison with the contour p lot of the 12.5 //m deconvolved image obtained at Keck 2 (fSock et al. 2000i) and the 
5GHz radio map of iGallimore et al.l (120041) (grey scale). Green dashed lines correspond to the more extended component. The 
blue ellipse corresponds to the elongated one at RA. ~ -4°, of size (< 140) mas x 1187 mas. The orange innermost component 
represents t he inner uniform disk < 85 mas. This component is directly associated to the dusty torus resolved in interferometry 
with MIDI dJaffe et al.lbOOllPoncelet et al.ll2006l). 



Fig.|4|i. The knot sizes are roughly equal to the diffraction limit 
of the telescope so they must be understood as upper size limits. 
In our case thanks to the super-resolution brought about by the 
interferometric analysis technique, the elliptical disk is thinner 
Moreover, the deconvolution process of iGalliano et al.. (.2005,) 
may have emphasized the intensity fluctuations in the structure 
just enough to produce the knots. On this basis, our ellipse must 
be viewed as an average model of these knots. Therefore, it is 
difficult to say which of the two source representations is the 
best, but it should be noted that the superimposition of the el- 
lipse and contour plot from deconvolved 12.5 /^m images from 
Bock et al. (l2000l) is very stricking (see right panel of Fig. |4]i. 



Our model does allow us to interpret a fraction of the flux of 
the central knot as part of the elongated structure, whilest the 
remaining fraction belongs to the central core (left unresolved 
by the image). Flux ratios between the different components 
are discussed in Sect. 13.2.41 

Examining the other wavelengths provides strong evidence 
for a correlation between the NS extended core and HST [OIII] 
clouds A and B (the loca tion of HST cloud B is marked in 
Fig. g] ICapetti etal]|l997h located at (y.'\ S to N of the core 
respectively. The NIR observations of iGratadour et al.l (120061) 
an d fainter levels of flu x observed on the deconvolved images 
of lGalhano et al.l (|2005|) suggest no depletion of dust in the re- 
gion between the core and knot NEl. 

The silicate feature evolution from the core up to + 3" 
outwards has been i nvestigated bv | Rhee & Larkin (l2006l) and 
Mason et all (|2006|) . iRhee & LarkinI (l2006l) interpret this trend 



as dust emission from the inner edge of the dusty torus. If this 



is correct it would mean that the torus should extend up to 1" 
in the EW direction, however neither interferometry nor VISIR 
observations have shown this to be true. Indeed molecular dust 
has been found inside the io nization cone up to 600 mas toward 
the N of the central en gine ( Mason et al. 20061) . Thus drawing 
upon the conclu sion of Mason et alj (|2006|) and using the the- 
oretical work of iDopita et al.l (l2002b . we explain our elliptical 
structure and the correlation between the NIR, MIR, and UV 
emissions in the nucleus of NGC 1068 in terms of unresolved 
small clouds and of molecular dust irradiated and photoion- 
ized by strong X/UV radiation fields produced by the central 
engine. Dust at the front is directly heated by X-rays then pho- 
toevaporates and generates a hot plasma of dust. This in turn is 
radiatively accelerated away from the central engine and emits 
in both NIR and MIR. In addition an envelope of ionized gas 
is produced that would account for the narrow emission lines 
in the UV and optical ranges. It would also explain the dis- 
tribution of [OIII] clo uds. This interpreta tion is strengthened 
by HST spectroscopy dGroves et alj|2004l) that shows that the 
CIII], CIV, and Hell UV line profiles have a pronounced blue 
wing within this region. This is associated with the proper mo- 
tions of velocities up to 3000 km.s"', which are indicative of 
NLR clouds being pushed toward us. Fits of these line ratios 
point towards photoionization as the dominating mechanism 
in comparison with excitation by radio-jet/interstellar medium 
shocks. Thus, NLR dusty clouds observed in NIR, MIR, UV 
must be at the front of the ionization cone and accelerated to- 
wards us in a radiative wind (see Fig. |5]i. The thinness of the 
ellipse can then be explained either by the ionization cone hav- 
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NIR extension (Gratadour et al £006) 
propagation of the shook wave ? 



Radio let (Gallimore et al. £004) 



Dusty molecular olouds (gradient of yellcm/ 
to red), individually unresolved in the IR + 
tail of dusty plasma (In blue) responsible for 
emission lines: [Olll] (Capetti et al. 1397), 
Clll], CIV, Hell (blue wing on line profiles 
indioative of velodty up to 3000 km/s from 
the systemic velocity; Groves et al. £004) 



Elliptical structure (<140) mas x 1187 mas 
traced by the speckle analysis of VISIR 
burst mode Images 




Plasma of dust radiatively accelerated (up 
to 3000 km/s from the systemic velocity; 
Groves et al. 2004) + dusty [Olll] tail = 
photoionlzation model from Dopita et al. 
(2002) 

+ IR observations from Galliano et al 
(2005) and Gratadour et al. (2006) 



Front shock with a dense cloud of the NLR: 
bending of the radio jet from P A =-4° to 
P.A. = 31 ° (Gallimore et al, £004) 



Clumpy molecular dusty torus, typical sizes 
from Jaffe et al. (£004) and Poncelet et al. 
(£006) + onentatlon {'■ 46°) from recent 
MIDI observations (Jaffe et al, £006) 



H20 maser disk (Green hi II &Gwinn 1397) 



Fig. 5. Picture summarizing the multi-wavelength structures observed toward the north in the nucleus of NGC 1068, as dis- 
cussed in Sect. 13.2.21 Amo ng ot h ers, th e representation of the radio jet (which bends from PA - -4° to P.A. = 31°) cornes from 
observations of Gallimore et al.l (|2004. H2O masers from iGreenhill & GwinnI (119971) . the dusty torus from iJaffe et al.l ( 120041) . 
Poncelet et"al1 ^200^ andljaffe et all (l2006l) . MI R clouds from the p resent study and fro mlGalliano et al.l (l2()03b . NIR emission 
from lGratadour et all (|2006|) . [QIIIl clouds from lCapetti et al.l ( Il997l) . and UV lines from lOroves et alj (|2004 . The yellow color 
was chosen for heated dust, the gradient of yellow to red color is used for decreasing dust temperature (typically ranging between 
300-500 K) and blue color for [OIII] UV clouds, for which Hell, CIII], and CIV emission-line profiles have a prononced blue 
wing indicative of proper motion up to -3000 km.s"', meaning that they are moving toward us. The blue elliptical structure 
derived from the analysis of VISIR BURST mode images is discussed in Sect. 13.2.21 



ing a very small opening angle < 50° (i.e. the angle between the 
two frontiers of the ionization cone), meaning that the open- 
ing angle of the torus would be > 150° (i.e. the angle between 
nothern and southern walls of the torus) or by the motion of the 
small unresolved dusty clouds being frozen in that of the radio 
jet. In both cases, their propagation would have to be confined 
to a thin region, such as what is observed. 

The effe cts of jet-NLR cloud s interactions are weak on the 
UV spectra dGroves et al.l 120041) so could not account for the 
high color temperature observed in the NIR ra nge so far away 
from the central engine dGratadour et al. I I2OO6I) . However, it is 
stricking that the elongated ellipse is aligned with the radio jet 
close to the core. There is actually some evidence for local jet- 
NLR interactions, such as the bending of the je t axis initially at 
PA = -4° to PA = 31° dGalhmore et alj|2004b . This would in- 
dicate a direct collision betwe en the radio jet and a dense NLR 
cloud- labeled HST-C in UV dCapetti et al.ll997l). NEl in MIR 



dGalliano et alj|2005h . and IR-Ib in the NIR dGratadour et al. 
200d) A twisting of structures involving t hese two PA is found 



in the deconvolved M IR and NIR images dGalliano et al ]|2005l 



Gratadour et al .l200d) . The peculiar distribution of knots above 



gation in the dusty pl asma medium as already discussed by 
Gratadour et al.l d2006l) . The global picture of this scenario is 
given in Fig.|5] 



3.2.3. More extended structures 

The extended component of the model of the VISIR visibili- 
ties has been used to account for the most extended contribu- 



cloud NEl may be the signature of the shock-wave propa- 



tionsto the flux. Given the deconvolved image of lGalliano et al 
d2005b . it includes knots NE2, NE3, NE4 and SW2 (see the 
middle panel of Fig.lH. The geometry we considered is clearly 
too simple, but a detailed description of the extended environ- 
ment of the core of NGC 1068 goes beyond the scope of the 
present study. The important information to consider here is 
the flux ratios between the total flux from the extended envi- 
ronment and from the inner structures. 



3.2.4. Flux ratios between components 

We compare the flux ratios values from our global modeling, 
equal to 3 and 1.4 between the different components, with 
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observations of Bock et al. I (I2OOOI) and iGalliano et all (l2005b . 
First, our inner component is associated to the source of flux 
labeled b in Fig. 5 of .Bock et al. (.2000) and the ellipse (re- 
ferred to elongated component in the following) to those la- 
beled a and c. According to their mean fluxes, we deduce 
^coi-eZ-Feiongated = 2.92. This ratio is in very good agreement 
with our best fit. 



Then, using Table 1 of IGalliano et al.l (l2005h and associ- 
ating both our elongated component with the SWl and NEl 
knots and our extended uniform disk with the remaining SW 
and NE knots, we have computed the following flux ratios: 

~ ^core /^elongated ~ 2.02 

~ ^elongated /^extended ~ 1.72. 

Given the uncertainties on photometry measurements in de- 
convolved images and the si mplicity of our model , and consid- 
ering that the structures of IGalliano et al.l ( l2005b are large in 
size compared to the diffraction limit of the telescope, these ra- 
tios are in good agreement with those derived from our study. 
Since the elongated structure in the EW direction and the cen- 
tral core are not resolved by the telescope, both fluxes are bi- 
ased, and this can explain part of the discrepancy between the 
two estimates. 

We conclude that our simple interferometric model has cap- 
tured most of both the sizes and relative intensities of the object 
structures when compared to the deconvolved images. It offers 
extra information thanks to the higher resolution achieved, both 
on the core associated to the dusty torus and on the width of the 
NS elongated structure. 

4. Comparison with M\D\ observations 

4.1. How to compare VISIR and MIDI visibilities ? 

One of the important aims of the study is the link between low 
spatial-frequency visibilities derived from VISIR images, with 
long baselines dat a points obtained with MIDI and the radiative 
transfer model by iPoncelet et al. I (l2006h . Visibilities extracted 
in part |2| correspond to a field of view for VISIR of several 
arcseconds, whilest the field of view for MIDI is limited to a 
slit of (y.'6x 2" oriented at - 30°. In addition, the emerging flux 
is mainly limited to what comes from the core by a ~ 0'.'6x 0'.'6 
mask oriented at - 30° centered on the nucleus (left panel of 
Fig.lHi. 

In processing images such as the conditions of MIDI obser- 
vations we first unbias spectra from the pixelization of VISIR 
images, and secondly we take the mask of MIDI into account. 
To correctly position the mask on the flux maximum, we have 
to keep the spatial information of images contained in the phase 
of spectra (in addition of the unbiased moduli). New spatial 
spectra of NGC 1068 and the calibrator are then calculated via 
the following equations: 



5nGC1068(«,v) = |VnGC1068 (m,v)| X e-'fNCCK 
5calib(",v)= |ycalib(",v)| X e"' '^'^'^ 



(5) 
(6) 



where IVNGCiOds (",v)| and IVcaiib (m,v)|, Ongcio68 and Ocaiib 
are respectively the moduli unbiased from noise (as calculated 



in Sect. |2]i and the phases of the spatial spectra of NGC 1068 
and of the calibrator star 

The bias produced from pixelization effects is corrected by 
dividing moduli by the pixel transfer function in the Fourier 
plane, the pixel being considered as a simple filter with a width 
of 75 mas (i.e. the size of one pixel) in the image plane. For 
oversampling, original 32 x 32 maps are padded out with ze- 
ros after the cut-off frequency of the UT (i.e. D/A), to achieve 
256 X 256 maps. The resulting maps produce new pixels an- 
gular size of ~ 9.4 mas.pixels ' . The MIDI mask coded onto a 
256 X 256 image is Fourier transformed. Spectra of NGC 1068 
and of the calibrator star are convolved with the spectrum of 
the mask, written as 



5mask(",v) = IVtnask («, v)| X e 



-i Dm 



(7) 



Low spatial-frequency visibilities derived from VISIR im- 
ages as obtained in the same conditions of MIDI data acquiere- 
ment are given by 



\V(u,v)\ 



\S 



NGC 1068 



(m, V) * 5tnask(M, V)| 



l'Scalib(M, V) * 5mask(M, V)| 



(8) 



where * stands for the convolution. 

MIDI observations made in 2003 correspond to baselines 
oriented at 0° and -45°. Therefore only cuts along these two 
orientations of low spatial-frequency visibility maps have been 
retained. For individual exposures, signal-to-noise ratio is low 
and phases are noisy. Therefore error bars on visibilities are not 
computed from the statistical distribution of individual mea- 
surements as in Sect. |2] but from the reduced of fits at 0° 
(see Fig.|7]i, which has been forced to be equal to 1. It is equiv- 
alent to considering the model as correct and to interpreting 
discrepancies with the data as errors. This leads to a qualitative 
estimate of the error of ~ 0.026, which has been applied to all 
visibilities between and 8 m. The comparison with MIDI data 
points is presented in Fig.|6l 

4.2. The link between MIDI and VISIR observations 

Figure |6] is the comparison between a) low spatial-frequency 
visibilities from VISIR images considering the MIDI mask, b) 
high spatial fr e quenc y ones from MIDI, and c) the model of 



Poncelet et al.l (120061) that reproduces MIDI data. The unex- 



pected sharp fall of low spatial frequencies visibilities demon- 
strated how constraining these measurements are for modeUng 
long baseline interferometric data. 

On the one hand, short baselines visibilities are well- 
reproduced by a uniform disk model of sizes ~ 450 mas and 
~ 350 mas along directions at 0° and -45°, respectively. On the 
other hand, high spatial-frequency visibilities from MIDI are 
represented by a radiative transfer model insid e a dusty layer 



ranging in size between 35 mas and 82 mas ( IPoncelet et al 
I2OO6). The initial attempt to fit both sets of data simultane- 
ously was associated with the model of the dusty layer of 



Poncelet et al.l (120061) being surrounded by a uniform disk. The 



only free parameters were the size of the uniform disk and the 
flux ratio between the disk and the dusty layer. This model 
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Fig. 6. Comparison between high spatial-frequency visibihties 
from MIDI at 12.8 fim (data points at 45 and 78 m) and low 
spatial frequency ones derived from VISIR BURST mode im- 
ages including the mask discussed in Sect. |4] (points between 
and 8 m of baseline). The model of the two uniform disks 
{black and grey solid l ines) has to be c ompar ed with the radia- 
tive transfer model of Poncelet et al. ( 20061) applied to MIDI 
data {black dotted line). A third component has been added to 
the one already revealed by the radiative transfer, to partially 
account for the steep fall of the visibility at short baseline. Size 
of this component is ~ 600 mas and ~ 300 mas along directions 
oriented at 0° and -45° respectively. 



agrees with the strong fall of visibility at low projected base- 
lines but fails to account for the faint level of visibility at high 
projected baselines (see Fig.|6ll. This forces us to attempt a new 
fit in w hich we have to cons ider the parameters of the dusty 
layer of iPoncelet et al. I (l2006l) as free. This is not yet possible 
according to the limited amount of BURST mode data sets that 
are currently available, as the wider spectral range that is re- 
quired. 

Nevertheless, to be able to quantify the fraction of extended 
emission entering the field of view of MIDI, compared to the 
emission of the compact core, we fit both low and high spatial- 
frequency visibilities together in both directions using a simple 
model of two uniform disks, in which free parameters are the 
sizes and the flux ratios between the two disks. The minimum 
value of the mduc&d is 3.42. This corresponds to an unre- 
solved inner disk of upper size limit ~ 21 mas in each direction. 
The extended component has a size of 479 + 17 mas along the 
NS direction and ~ 350 + 9 mas along PA = -45°. The flux ra- 
tios between the two components are found to be 12 + 5 and 
6 + 2 in these two directions (see fits in Fig.|7]i. 



4.3. Discussion 

This study demonstrates that there are two main structures 
appearing in the MIDI field of view, a compact one of size 
~ 20 mas surrounded by an extended and slightly elongated 
one (the size ratio between -45° and 0° is ~ 27% for this sec- 
ond component). The first structure is associated to the dusty 



Fig. 7. Comparison between high spatial-frequency visibilities 
from MIDI at 12.8 /im and low spatial frequency ones derived 
from VISIR BURS Tmode images including the MIDI mask as 
discussed in Sect. |4] Solid black and dashed grey lines corre- 
spond to a model of two circularly uniform disks applied along 
both directions at 0° and -45° independently. Optimum parame- 
ters are: 0i(O°) and 0i(-45°) < 21 mas, 02(0°) = 479 +17 mas, 
02(-4°) = 450 + 9mas, 77(0 deg) = 12 ± 5 and 77(-45°) = 6 + 2. 



layer highlighted so lely by the MIDI data (iJaffe et al.1 12004 



Poncelet et al. 20061) . Sizes are slightly different due to the 



fact th at it is a more simp l istic rn odel than those of iJaffe et al 



( |2004 and lPoncelet et al.l ( I2OO6I) . In addition, the direct envi 
ronment surrounding the compact core (of size ~ 500 mas) that 
enters the field of view of MIDI was not taken into account by 
previous studies. According to the flux ratio between compo- 
nents, it contributes to more than 83 % of the total flux seen by 
MIDI. 

Constraints o n temperatures of the dusty torus derived from 
only MIDI data ( Poncelet et al. 20061) are not biased by this 
effect. Indeed they come from an additional information be- 
ing the spectral energy distribution supplied by MIDI. They 
are found to be on the order of ~ 300 K, a value supported 
by the continuum slope of the central 0.4" N-band spectra 
jMason et al. I l2006h . That temperatures are so low so close to 
the central engi ne is strong evidence of t he clumpy nature of 
the dusty torus ( Ehtzur & Shlosmanll2006 ). 

The compact central component is already traced with the 
study of full field images presented in Sect. [3j it is directly as- 
sociated to the innermost component < 85 mas. The extended 
component, of sizes 479 + 17 mas and 350 + 9 mas along 
PA = 0° and -45°, is not entirely related to the elliptical NS 
elongated ellipse traced in Sect. 13.1.21 since this last one does 
not completely enter in the field of view of MIDI. However, 
according to the d iscussion in Sect. 13.2! and to the photoion- 
ization scenario of iDopita et al.l (120021) . this elongated compo- 
nent could be associated to small dusty clouds heated by the 
radiation field from the central engine and p hotoevaporating . 



The NIR and UV counterparts of these clouds jGratadour et al 



20061: ICapetti et al.lll997l;lGroves et all2004 and the evolution 
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of the 9.7 yum silic ate emission toward the north of the core 
dMason et al. I l2006h is evidence for the presence of dust in- 
side the ionization cone at distances less than 300 mas from 
the central engine (i.e. half of the MIDI mask). The present 
study supports this fact and demonstrates the need to consider 
the close environment of the dusty torus for modeling interfer- 
ometric data. 

The small set of MIDI and VISIR BURST mode data cur- 
rently available only allows the use of simple geometrical mod- 
els for comparing low and high spatial-frequency visibilities. It 
is strongly necessary to apply this approach at different wave- 
lengths to get a better description of the source of mid-IR emis- 
sion seen by MIDI, associated with the dusty torus, and linked 
with its immediate surroundings. 

5. Conclusions 

The speckle analysis of the 12.8 /im BURST mode VISIR 
images constrains the upper size limit of the MIR core of 
NGC 1068 to 85 mas. This is well below the diffraction limit 
of the UT which demonstrates the achievement of the BURST 
mode of VISIR. The NS elongated component is described here 
as an elliptical structure oriented at PAs = -4° (such as the base 
of the radio jet) and of siz e (< 140) mas x 1 18 7 mas. Following 
the scenario proposed bv lDopita et al.l (120021) . it is interpreted 
as small dusty clouds individually unresolved, distributed in- 
side the ionization cone, and photoevaporating. This explains 
the correlation of observations at different wavelengths. 

The use of VISIR BURST mode images for the compari- 
son with MIDI data reveals the strong and unexpected fall in 
visibility between and 6 m of baseline. The use of simple ge- 
ometrical models accounts for the two data sets simultaneously 
and distinguishes between two main components: a compact 
one ~ 20 mas asso ciated to the dusty torus traced by inte rfer- 
ometry with MIDI djaffe et al.ll2004l: IPoncelet et al]|2006l) and 
an extended one (< 500 mas) entering in the field of view of 
MIDI. This surrounding environment, partially related to the 
elongated elliptical structure, contributes to more than ~ 83% 
of the flux emitted by the core. It is therefore mandatory to take 
it into account when modeling the MIDI data. By doing this 
we do not challenge the low temperatures of the dusty torus 
constrained by the previous analysis of MIDI data. Besides, it 
reinforces the scenarios of clumpy torus. 

According to the small set of VISIR BURST mode and 
MIDI data available up to now, we reach the limits of sim- 
ple descriptions of sources of mid-IR emission in the AGN of 
NGC 1068. To further describe the complexity of the compact 
dusty core and the deeper link with its direct surroundings, the 
same processing of VISIR BURST mode images has to be per- 
formed at other wavelengths. This will be possible when the 
BURST imaging mode of VISIR is open for observations at the 
VLT. Moreover, as it shows how constraining low spatial fre- 
quency visibility points are, this study underlines the need to 
further fill the gap of interferometric data between 8 and 45 m 
of projected baseline. On the one hand, baselines between 30 
and 45 m are reachable with the VLTI and the UTs. On the 
other, since most of the ATs are now installed at Paranal, it will 
be interesting to attempt observations with these 1.8 m pupil 



telescopes to bridge the gap from 8 to 23 m to allow for more 
accurate multi-scale modeling of the inner part of NGC 1068. 
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